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We propose a Z2-stabilised singlino (x) as a dark matter candidate in extended and ^-parity violating versions 
of the supersymmetric standard model. % interacts with visible matter via a heavy messenger field 5, which 
results in a supersymmetric version of the Higgs portal interaction. The relic abundance of % can account for 
cold dark matter if the messenger mass satisfies Ms ,$ 10 4 GeV. Our model can be implemented in many realistic 
supersymmetric models such as the NMSSM and nMSSM. 
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I. INTRODUCTION 

It is well established that visible matter is not sufficient 
to account for the observed structure of the Universe. This 
implies the existence of non-baryonic dark matter (DM). 
Global fits of cosmological parameters can accurately deter- 
mine the density of DM, albeit indirectly. Measurements of 
the anisotropy of the cosmic microwave background radiation 
(CMBR) and of the spatial distribution of galaxies give for the 
density of DM {J 

£W* 2 = 0.106 ±0.008. (1) 

Identifying the nature of dark matter is a major goal in as- 
troparticle physics. Many particle physics candidates have 
been proposed in both supersymmetric (SUSY) and non- 
supersymmetric extensions of the standard model (SM). In ei- 
ther case the stability of DM is ensured by imposing a global 
symmetry. The simplest global symmetries considered are Z2 
and t/(l); see for example (j, l2l U, |5l |6j |7[] . 

In low energy effective SUSY theories the symmetry is usu- 
ally /^-parity, (— \^( 3B+L+2S \ which is imposed to conserve 
baryon (B) and lepton (L) numbers. As a result the stability 
of proton is ensured. It turns out that R is +1 for all SM fields 
and -1 for their superpartners. Thus /^-parity, which is a Z2 
symmetry, protects the decay of lightest SUSY particle (LSP) 
to SM particles. As a result the LSP is a good candidate for 
DM within minimal SUSY standard model (MSSM) and its 
extensions as long as the conservation of 7?-parity is ensured. 

However, B and L are accidental global symmetries of SM. 
Thus it is not clear a priori that B and L are conserved within 
the MSSM. If B and L are violated then /^-parity is not con- 
served. Non-conservation of /^-parity is one way to generate 
small neutrino masses [8], which provide solid evidence for 
phyiscs beyond the SM. Moreover, if /^-parity is violated then 
leptogenesis is possible [9], which explains the small mat- 
ter anti-matter asymmetry (o(10~ 10 )) required for successful 
Big-Bang nucleosynthesis. However, within the MSSM and 
its extensions there is no well-motivated particle physics can- 
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didate for DM in the presence of /^-parity violation . 

In the following we will explore an alternative possibil- 
ity for the DM candidate in SUSY models, irrespective of 
whether /^-parity is violated or conserved, by introducing a 
new Z2 symmetry and additional singlet fields. Singlet exten- 
sions of the MSSM are often considered to ensure that the /j 
parameter is at the electroweak scale 111 111 . The prime among 
them are the NMSSM (the Next-to-Minimal SUSY Standard 
Model) and the nMSSM (the nearly-Minimal SUSY Standard 
Model). In such models, if /^-parity is conserved then the DM 
candidate can be an R -parity odd singlino lfl2tl . Here we pro- 
pose an alternative SUSY DM candidate: a Z2-odd singlino 
© which is stable without requiring /^-parity 2 . 

Beyond considerations of /^-parity, Z2-singlino dark mat- 
ter is interesting as a SUSY implementation of gauge singlet 
dark matter. Gauge singlet scalar dark matter interacting via 
the Higgs portal 11411 was first discussed in detail in Jj], with 
a further study presented in and an earlier analysis given 
in J6J]. With the advent of the LHC, Higgs portal couplings 
to hidden sector particles are of considerable topical interest. 
The superpotential coupling we will consider here is the nat- 
ural extension to SUSY of the Higgs portal concept. How- 
ever, it is necessarily non-renormalisable are a consequence 
of SUSY, pointing to the existence of further new particles at 
the TeV scale. 



II. MODEL FOR Z 2 -SINGLINO DARK MATTER 

A. R-parity conserving SUSY 

We extend the MSSM by adding a chiral superfield % and a 
messenger field S. We also impose an additional Z2 symmetry 
under which % is odd, while all other fields are even. The full 



In supergravity (SUGRA) theories, the gravitino can account for dark mat- 
ter in certain regions of parameter space since its coupling with matter 
fields is suppressed by the Planck scale llOll . 
2 A different Zi -singlino dark matter model, which is based on a broken U ( 1 ) 
gauge group, was presented in fiUl . 
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superpotential is 

W = Wmssm + hS%% + X 2 SH u H d - 

z, z, 

where 

Wmssm = h^L^Hd + hfjQ iU c H u + hfjQid c H d +fjH u H d . (3) 

In this case the effective superpotential after integrating out S 
becomes 



M- 



X „2 



W = Wmssm + 



fX 2 H u H cl 
M S 



(4) 



where f = X\hi- The term with coupling / is the natural gen- 
eralisation to SUSY of the Higgs portal-type coupling to % 
scalars of the form X '%H^H IU4ll . However, SUSY implies 
that the Higgs portal interaction is now non-renormalisable. 
The Lagrangian terms involving the interaction of % scalars 
and fermions, to order 1 /M$, are then 



C. Gauge singlet dark matter 

Both the scalar and fermion components of the x superfield 
are stable due to the Z 2 symmetry and therefore the lightest 
of these will be a potential DM candidate. In most cases 
the lightest component will be the fermion, the Z2-singlino 
% since the scalar component will gain additional mass from 
SUSY breaking. We will therefore focus on the Z2-singlino 
as the DM candidate 3 . Its relic abundance will then be deter- 
mined by the following scattering processes: 



XX 
XX 
XX 



MSSM fields 

xh 

MSSM fields . 



(9) 



The latter two processes will be negligible due to Boltzmann 
suppression if the x mass is large compared with the % mass. 
We will assume this to be the case in the following. There- 
fore we will only consider the first class of processes when 
calculating the relic abundance of %. 



~H D \ M t\ 2 xh + M t x-x 



-r^XX^HMd 

Ms 



~H u H d X.% + h.c. 
M s 



■0(l/M 2 s ), 



(5) 



where x denotes the scalar and x the two-component fermion. 



B. fl-parity violating SUSY 

The superpotential involving /^-parity non-conserving inter- 
actions is: 



M 7 , X LiH u 



(6) 



where 



W fp = hjkULjll + X ijk LiQjd L k + X i}k u c id c jd c k + lhUH u (7) 

is the R -parity non-conserving superpotential in MSSM. The 
R -parity violating terms in the Lagrangian involving the inter- 
action of x, to order 1 /Ms, are then given by 



H 3 \M x \ 2 xh+M x X-X 



hi My + ~ 

- — ^XX LiH u 



Ws^ 1 '-"" + Ws lHlllLi + Ws 11 ' 1 -" 11 



hi - 

JT L i H uX-X + h.c. 
M s 



where L, is the slepton doublet. 



■0(1 /Ml), 



(8) 



III. RELIC ABUNDANCE OF Z 2 -SINGLINOS 

In this section we calculate the relic abundance of %. We 
first calculate the scattering cross-section times relative veloc- 
ity for annihilation processes to MSSM final states. 

After electroweak symmetry breaking there are five phys- 
ical Higgs scalar degrees of freedom. In this letter we will 
consider the physical Higgs scalars to correspond to gauge 
eigenstates when calculating the cross-sections, with all Gold- 
stone bosons coming from H d . The physical Higgs scalars are 
assumed to have a common mass Mh- In addition, we will 
consider the gaugino and Higgsino gauge eigenstates to cor- 
respond to mass eigenstates with a common neutralino mass. 
A more general analysis will be presented in future work. 

H , H„ 




\H d ,L, 

FIG. 1 : The four point annihilation of %% t0 Higgs an d sleptons in 
MSSM 



In the non-relativistic limit the contribution to the total an- 
nihilation cross-section times relative velocity of XX annihila- 



' There may be regions of parameter space where the SUSY mass M z is 
close to the SUSY breaking mass terms, in which case the scalar % could 
be the lightest component. We will return to this case in future work. 
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tion to Higgs and sleptons (Fig. 1) is given by: 



(12) 



i M- 

<<nM) = — ^(l+v4/2) 



4%sMj 



Mr 



(10) 



The contribution of %% to SM fermions through /^-parity con- 
3>„ 






<I>a 




w,z 



FIG. 3: Mutual annihilation of XX to gauge bosons through Higgs 
and sleptons. Here * a = H®,H®. 



The contribution of %% to W -bosons (Fig. 3) is given by 



FIG. 2: Mutual annihilation of x to SM fermions through Higgs and 
slepton. Here If, F=Q h L u H d and F c =u C j,7j,Tj. 



serving interactions (Fig. 2) is given by 
i / M-\ 

<f»\**\) = ^ v (1+v - /2) 



Wd) 

M S 



2\2 



(l-^) 2 



M S J ' (, Ml- 2 



2M± 



u |2 



(>-?)' 



,7(g»)V|y |2 ll ~~ 



2H 



(l-*y 



(11) 



The contribution of %% to SM fields through R -parity violating 
interactions (Fig. 2) is given by 

i / M-\ 

«*M> = (1+^/2) 

hi(H u )\ 2 lue A s s 



M s 



M- 



2 \ 2 



+ 



h,(H u ) , 2 



1 



1-2?)" 



1 - 



M7 



+ 



1 __^L 



(-.M>-^(. + ^)(2 + ^£ 



+ 



(mi 



2 \2-| 



2\2 



(13) 



where v = ((#„) 2 + {H d ) 2 ) l l 2 = 246 GeV, while the contribu- 
tion of to Z-bosons (Fig. 3) is given by 



hM .'5 (lt4 „( !t a: 



•2\ 2 



^) V) 2 (^) 



(^) Z (^) 2 ( ? 

(,-M 2 ) 2 



+ 



(14) 



The contribution of XX to sparticles and MSSM Higgs bosons 




B 



FIG. 4: Mutual annihilation of XX t0 sparticles and Higgs. Here 
*a = H^Hj,!^ and A, B stands for the sparticles and Higgs. 
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(Fig.4) is given by 



(<*6 V re l 



J_Mf 
471 s 



(l+V r 2 el/2) 



ff(H d )Y , . . 



( Wu) 

(s — M~ ) 2 fjj 
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.? 






_M| 


s 


s 



(15) 



where is the mass dimension coupling at the tri-linear 
scalar vertex. Finally, the contribution of to gaugino and 





FIG. 5: Annihilation of XX t0 gauginos and fermions. Here <I> a 
H°,Hj,L1, F = K,m,Li and X = W,B. 



fermion (Fig. 5) is given by: 



IV. CONSTRAINTS ON ^-PARITY VIOLATING 
INTERACTIONS AND IMPLICATIONS FOR % 
ANNIHILATION 

Before estimating the relic abundance of % let us briefly 
discuss the constriants on /^-parity violating interactions (O 
in the MSSM In MSSM there are three types of tri- 

linear ^-parity violating couplings: A,,-^, X j j k and While 

Xijk is antisymmetric with respect to i and j, X"- k is anti- 
symmetric with respect to j and k. Thus the /^-parity vio- 
lating interactions in general add 45 extra parameters to the 
MSSM. These couplings are severly constrained by the non- 
observation of certain physical phenomena. In particular, the 
product XX < 10~ 9 comes from the stability of proton. Simi- 
larly, non-observation of n — h oscillations gives the constraint 
X" < 1(T 5 for m = 100 GeV, where m is the SUSY breaking 
mass. The X and X couplings induce a Majorana mass for 
three generations of neutrinos. The electron neutrino mass 
then gives the constraint X.X < 10~ 3 for m = 100 GeV. Neu- 
trinoless double beta decay gives the constraint X < 10~ 4 . 
Thus we see that these trilinear couplings are necessarily small 
in comparison to /^-parity conserving couplings in the MSSM. 
Therefore, the annihilation channels of through these trilin- 
ear R -parity violating couplings are necessarily small in com- 
parison to the /^-parity conserving couplings. 

There is a bilinear term PjLiH u in the R -parity breaking su- 
perpotential. However, one can show that by making a SU (4) 
rotation /U ( L,H M can be rotated away Illql . leaving only the bi- 
linear term iiH u li ( i which is /^-parity conserving. Therefore 
the presence of such a bilinear term in the /^-parity breaking 
superpotential does not contribute to any extra annihilations 
of%- 

In what follows we neglect all annihilation channels of jx 
to MSSM fields involving /^-parity violating couplings X, X 
and X . However, we note that the new R -parity violating cou- 
plings hj are not necessarily small. When estimating the relic 
abundance of % we will consider only those /^-parity violating 
channels involving the couplings hi. 



M s 



I _ H u X 



- A — * 



( f " ' ' * 

\ M S 



M\ Ml 
9 I 1 - 1 L 

M s ) I u. 



(16) 



V. DENSITY OF Z 2 -SINGLINO DARK MATTER 

The relic abundance of % can be calculated by solving the 
Boltzmann equation: 



dn 



—^ + 3njH = - (o ann | v rd | ) - n| c 



eq2 



(17) 



where (c ann \ v re i ) is the thermal average of the %X annihilation 
cross-section times relative velocity, with a a ,m — T^i^iJ = 
1 , 7, and nj is the number density of %. The equilibrium den- 
sity of non-relativistic % particles is 



eq ~ 

'¥ = 2 



MjT 
2n 



-i 3/2 



-Myr/T 



(18) 
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With / = n^/T 3 , Eq. (Qj) becomes 



df 
dT 



K 



(19) 



where = «l q /7 3 and - 



[4n 3 g(T)/45M 2 pl ] 1/2 . The dcn- 



sity can then be calculated using the Lee-Weinberg approxi- 
mation 11711 . The freeze-out temperature, 7b, is defined by 



dfeq 

~dT 



in | Vrel | ) ,2 

K Teq 



(20) 



To obtain the present density Eq. ( fT9l ) is solved from 7/j to the 
present with f eq = on the right-hand side and with /(7b) = 
feq{To)- The freeze-out temperature can be described by a 
dimensionless parameter zd ~Mj/Td- Solving Eq. (f20b gives 
for zd, 



Ad- 



ZD 



X 

T D 



In 



0.076- 



1 M f Mpi{o ann |y ref |) 
1/2 1/2 / _3_ 



(21) 



where = g(7p) is the effective number of relativsitic de- 
grees of freedom at Td- This implies that zd ~ 25. Solv- 
ing Eq. ( fT9] l with / eg = on the right-hand side and with 
/{Td) = feq (Td) then gives the number density at a lower tem- 
perature, 



nj{T) 



g(T) x 1.67 g y 2 7 3 z Q-^) 
A%M rt (a awi |v re/ |) ^ _ 



; 7<7 D 



(22) 

where we have included a correction for the change in the ef- 
fective number of relativistic degreees of freedom. Therefore 
the present contribution of % to the critical density of the uni- 
verse is 



Q.yh z w 1.1 x 10 y GeV~ 



ZD 



1 /2 

gj Mpi(o 

aim 



(23) 



where zd 1 is assumed. 

In the following we consider in the limits (i) s < 
A/|,Af|, and (ii) 5 > A/|, 
relativistic limit. 



M|, where 5 



4Mj in the non- 



(i) Small A%: s<M%,M? 



To focus on a definite example we set Mh = Mi = 150 
GeV and tanp = (H u )/(H<i) = 1 in the cross-sections. We 
assume that the mass of the other sparticles is 100 GeV. 
Since we assume that s < Mff,M^ in this case only <T2, 03 
and 07 will contribute to the relic abundance of %. We put 
/ = hi = 1; the results for smaller values can be obtained 
by rescaling Ms. In this case the allowed region in the plane 
of M s versus M% for £ljh 2 = 0.106 ± 0.008 is shown in Fig. 
©. It can be seen that for 15 GeV < My < 50 GeV, M s is in 



the range 1-3 TeV. The behaviour can be understood as fol 
lows. In the limit s < Mjr^M^ , the annihilation cross-sectioi 
Oann = 02 + 03 + 07 times relative velocity is of the form: 



\ v rel\) ' 



M~ 

M\ 



(24) 



where C is a dimensionful constant involving the VEV of H u 
and H c /. Therefore, smaller values of Mj require small values 
of Ms in order to keep £ljh 2 constant. 
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FIG. 6: Contour of %/j 2 = 0.106 ±0.008 is shown in the plane of 
My versus Mj. We have taken f = hj= 1 . 



(ii) Large Mj: s > M%,Mi 



We next consider s > M^,M?. We show the allowed re- 
gion in the plane of Ms versus Mj, corresponding to £ljh 2 
0.106 : 0.008. in Fig. ©. From Fig. © it can be seen 
that for Mj > 200 GeV, Ms is almost constant at around 
10 3 - 84 GeV = 6.9 TeV. This can be understood as follows. 
In the limit s > Mjj,M^ , the annihilation cross-section a ann = 

a, (i = 1 — 7) times relative velocity is of the form: 



\ v rel\) 



1 

Mj 



-CI 



1 



M \ M l, 



(25) 



For M Y > 200 GeV, the 



where C is a dimensionful constant 
effecive annihilation cross-section is dominated by the first 
term. As a result we get a constant value Ms ~ 6.9 TeV. For 
Mj < 200 GeV, the second term in the above equation dom- 
inates. In this regime, larger Ms is required to keep £ljh 2 
constant as My_ decreases, with a Higgs pole at M% = 75 GeV 
allowing much larger values of Ms over a small range of Mj. 
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In general, smaller values of Ms are possible by reducing / 
and hi, so the values shown in the figures should be considered 
as upper bounds on Ms, corresponding to large / and hi. 
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FIG. 7: Allowed region of %/i 2 = 0.106 ±0.008 is shown in the 
plane of Ms versus M%. We have taken f = hj = 1 . 



VI. CONCLUSIONS AND OUTLOOK 

We have discussed the possibility of Z2-singlino dark mat- 
ter in extensions of the MSSM. The dark matter communi- 
cates with visible matter through a heavy messenger field, S. 
As a result the interaction is suppressed by the mass scale Ms- 
For Ms £ 10 4 GeV the Z2-singlino can be cold dark matter 
for a wide range of mass, 15GeV < Mj < ITeV. (Larger val- 
ues of Ms are possible near a Higgs pole.) The possibility 
of dark matter in this case does not rely on the conservation 
of /^-parity. Thus the model is particularly important for the 
MSSM and its extensions, such as the NMSSM and nMSSM, 
when R -parity is violated. Non-conservation of R -parity is of- 
ten considered to give small neutrino masses, as required by 
the oscillation data, and for leptogenesis, a robust mechanism 
for the matter anti-matter asymmetry of the Universe. 

In the case of non-SUSY gauge singlet scalars interacting 
via the Higgs portal, direct and indirect detection rates are 
comparable with conventional weakly interacting dark matter 
candidates fill!. In the Z2-singlino case the coupling to the 
Higgs has an additional suppression factor as v/Ms, where v is 



a Higgs expectation value. Therefore we would expect signif- 
icant detection rates for Ms "Z 1 TeV. In this case the effective 
theory based on integrating out the S fields may not be appro- 
priate. We will return to the question of Z2-singlino detection 
in future work. 

The Z2 symmetry responsible for dark matter in this model 
can be a surviving symmetry (a discrete gauge symmetry) of 
a gauged U(\)' extension of MSSM. Such models are nat- 
ural in top-down scenarios when E(6) grand unified theory 
is broken down to the MSSM. A gauge origin of the Z2 is 
favoured by arguments which suggest that global symmetries, 
both continuous and discrete, are broken by non-perturbative 
gravitational effects [19]. In this case 7?-parity may be broken 
while a Z2 discrete gauge symmetry may account for SUSY 
dark matter. 

We have focused on the case of Z2-singlino dark matter pro- 
duced by conventional freeze-out from thermal equilibrium. 
There is, however, another possibility. In the case of non- 
SUSY gauge singlet scalar dark matter, when the mass of the 
scalar is entirely generated by the Higgs expectation value, the 
correct relic density is produced via decay of thermal back- 
ground Higgs bosons when the mass of singlet scalars is in the 
range 1-10 MeV [20]. This is the ideal range 1I2III22I1 for very 
long-lived dark matter particles to account for the 511 keV 
line obeserved by INTEGRAL |23]. In the Zz-singlino model, 
the singlino mass will be entirely generated by the Higgs ex- 
pectation value in the limit M z — > 0. We will consider the light 
Z2-singlino in a forthcoming paper 

Although we have considered dark matter particles which 
are Standard Model singlets, the model can easily be gener- 
alised, for example, to a SUSY version of the inert doublet 
dark matter model 0]. In addition, the messenger mass in 
the model can be greater than 10 4 GeV, in particular for the 
case where the Z2-singlino mass is close to a Higgs pole. This 
may allow the messengers to be associated with the messenger 
fields of a gauge mediated SUSY breaking model. 

The model we have presented here may be regarded as a 
SUSY generalisation of the Higgs portal concept. As such, 
we can expect the model to arise in the low energy effective 
theory of a wide range of SUSY particle physics models. 
Comment: While this pap er was in preparation a similar 
model was presented in 1I24I1 . 
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